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Abstract 

We propose and investigate a model for the coexistence of Superconductiv- 
ity (SC) and Ferroelectricity (FE) based on the dynamical symmetries su{2) 
for the (pseudo-spin) SC sector , /i(4) for the (displaced oscillator) FE sector 
and su{2) /i(4) for the composite system. We assume a minimal symmetry- 
allowed coupling ,and simplify the hamiltonian using a double mean-field ap- 
proximation (DMFA). A variational coherent-state (VCS) trial wave-function 
is used for the ground state : the energy, and the relevant order parameters 
for SC and FE are obtained . For positive sign of the SC-FE coupling co- 
efficient, a non-zero value of either order parameter can suppress the other 
one (FE polarization suppresses SC , and vice versa). This gives some sup- 
port to "Matthias' Conjecture" [1964] , that SC and FE tend to be mutually 
exclusive. For such a Ferroelectric Superconductor we predict: a) the SC 
gap A (and Tc) will increase with increasing applied pressure when pressure 



quenches FE, as in many ferroelectrics and b) the FE polarization |P| will 
increase with increasing applied magnetic field up to He, which is equivalent 
to the prediction of a new type of Magneto-Electric Effect in a coexistent 
SC-FE material. Some discussion will be given of possible relation of these 
results to the cuprate superconductors . 

Keywords: Superconductors , Ferroelectrics, Coexistence, Dynamical Sym- 
metry. 
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I.Introduction, and Background. 

This paper concerns an investigation of coexistence of Superconductivity (SC) and lattice 
Ferroelectricity (FE) based on a model which expresses the Dynamical Symmetry underlying 
the physics. The model combines the dynamical symmetry of the SC and FE sub-systems 
into that for the composite system. 

In Section II we review the algebra of the "s-wave" pairing BCS model of a super- 
conductor . It is well known that a dynamical symmetry su{2) algebra can be identified 
based on time-reversed electron pair operators . The mean-field reduced Hamiltonian will 
be an element in this algebra. By rotating the Hamiltonian in the space of the generators 
to a "diagonal " form, the energies , the eigenfunctions and then the expectation value of 
the SC order parameter in the ground coherent state, is obtained. 

In Section III we introduce a simplified algebraic model which has /i (4) dynamical 
symmetry for a displacive ferroelectric. It is a "displaced oscillator" model for the phonon 
soft-mode which has represening the soft transverse-optic (TO) mode. This Hamiltonian 
can also be transformed to "diagonal" form to give the energies, the eigenfunctions, and the 
FE order parameter in its ground coherent state. 

In Section IV we introduce the SC-FE coupling ,and we discuss interactions which 
will respect the gauge and inversion symmetries that will be broken. The Hamiltonian of the 
composite system including initially biquadratic interaction, is simplified using a "double 
mean field approximation ". The resulting bilinear Hamiltonian Hdmfa is in the direct 
product sm(2)®/i(4) algebra, and appears to be the simplest way that the two sub-symmetries 
can be joined. 

In Section V we show that although this Hamiltonian cannot be solved exactly, a 
variational solution can be found, by forming a trial eigenfunction analogous to the product 
coherent state. After carrying out the variational solution, the energy spectrum, and the 
eigenfunctions are obtained. 

In Section VI we calculate the expectation values of the SC and the FE order pa- 
rameters in the ground state of the coupled system. Our results for the order parameters 
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show that : the presence of one non-zero order parameter (e.g. spontaneous polarization) 
tends to suppress the other (e.g. superconductivity) and vice-versa in the case of positive 
sign (repulsive) of the coupling between the two subsystems. This leads to the prediction 
that the superconducting critical temperature can increase with pressure if FE is quenched 
by pressure , and is in agreement with experiment in sodium tungsten bronze. Another 
prediction is that the ferroelectric polarization will increase with applied magnetic field . 
This is a new type of Magnetoelectric Effect. These predictions are discussed in Section VI. 

In the final Section VII we discuss the results especially for the sodium tungsten 
bronze, and doped SrTiOs systems and we suggest the relation of the present work to the 
mechanism of high temperature superconductivity in the copper-oxide systems as previously 
proposed by Peter and Weger [Q based on the close proximity of a near- ferroelectric 
instability to the superconducting transition. 

We believe the model presented here offers a simple "generic" way to treat the two 
broken symmetries relevant to the problem: for SC, broken gauge symmetry ( the basic 
pair operators do not conserve number of electrons), for FE broken inversion symmetry 
(spontaneous TO phonon displacement is not i-invariant). In passing we will show later 
that , if certain couplings vanish , our model can reduce to other, well-studied models: the 
Jaynes-Cummings Model , and the Spin-Phonon model. Our model is more general than 
either of these two, and truncating our model to obtain either of them would result in losing 
relevant physics for our case. 

Apart from interest in new models for competing phases , there are several reasons 
for work on this problem at this time. As far back as 1964 , Matthias and coworkers 
reported superconductivity with Tc below 1°K in a sodium tungsten bronze NaxWOs with 
0.1 < X < 1. These authors remarked, the host crystals are isomorphic (at a; = 1 ) to barium 
titanate,so that it is "probable that they are also Ferroelectric, in the sense of developing a 
polar axis, similar again to to BaTiOs and WO3" . Matthias (1949,1967) and Matthias 
and Wood (1951) had confirmed the polar state for the sodium tungsten bronzes which 
they studied. Subsequent work of S.C.Abrahams et. al. [Q, and others, has also reported 
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on the structural phase transitions and the development of a polar axis in many of these 
compounds. The doped tungsten bronzes have also been studied as examples of the Mott 
and metal-insulator transitions 0. In undoped WO3 , five phase transitions have been 
identified: at 40, 65, 130, 220, and 260 °K 0, [0, and most recently such transitions were 
studied by by Aird et.al. Several recent reports of high temperature superconductivity 
Tc ~ 90° K in a sodium tungsten bronze NaoQ^WO^ system added considerable new 
stimulous for this work, since we can suppose that the system is in a polar state at the 
superconducting temperature , which allows ferroelectricity. 

It is worth recalling, too, that SC-FE coexistence, and competition played a role in 
motivating the work of Bednorz and Muller |]10| on the high temperature cuprate super- 



conductors. Even earlier, work on the "old " superconductors of the (3 — W structure like 
V^Si^Nh-^Sn^eic, where T^. ~ 23"/^, and a martensitic phase transition occurs in the same 
temperature range , gave rise to investigations on the possibility of a "ferroelectric metal" 
or a "polar metal" , and thus to the study of SC-FE coexistence ||Tl|, |T^, IT^, p5[] . 



A number of theoretical papers have already discussed microscopic models for the 
effect of lattice instability on superconductivity in the sodium tungsten bronze systems 
P^ , pS| . These papers have illuminated many aspects of the interplay between the 



structural deformations , such as rotation of underlying octahedaral units, and coupling with 
electron pairs. The present work looks at the same problem of coexistence of SC and FE 
from a dynamical algebra , or "Spectrum Generating Algebra [SGA]", point of view , which 
complements these detailed models . 

Study of SC-FE coexistence problems can be relevant to recent work by Weger 
and collaborators fl]], 0, on the mechanism of high temperature superconductivity in the 
cuprates. In that work the presence of a nearby FE instability close to the SC transition 
is related to the anomalously large ionic dielectric coefficient in the cuprates, which re- 
duces the electron-electron repulsion and can then lead to an enhanced net electron-electron 
attraction, producing higher T^.. 

The present work also relates to earlier dynamical symmetry investigations by Birman 



and Solomon 0, 0, |2[, [||], [||], l2|], for systems with multicritical 
behavior involving superconductivity and charge and spin density waves. Based on the 
mean-field models which had earlier been used to investigate superonductor-charge density 
wave , and superconductor- antiferromagnetic coexistence, these papers were the first to 
introduce models with S0{6) — 5*0(5) symmetry for Superconductor plus Charge Density 
Waves (SC-CDW), and Superconductor plus Antiferromagnetism (SC-AFM), as well as 
presenting and analysing a general SU{8) "Grand Unified Theory (GUT)" model unifying 
singlet and triplet Superconductivity, and Charge and Spin Density Wave cooperative effects 
||26|| . Recently , S0(5) and SU(4) , models for multicritical superconductor-antiferromagnetic 
behavior in the high temperature superconductors have also been studied by Zhang, Demler, 
Guidry , and others |2^, ^9 . 
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II. The su{2) Pseudospin Model For a Superconductor. 

The BCS theory for superconductors can be conveniently epitomized at a mean-field 



level by introducing the pseudo-spin su{2) algebra of the fermion pair operators ||30|, ^1 
p^ , [Q. Since this is well known, here we briefiy summarize the material needed for later 
reference. Some more notational and other details are given in Appendix 3. 

In the dynamical symmetry 5^(2) model for a superconductor we take the Hamilto- 
nian in a reduced mean field approximation as 

Hsc = Y.hk (1) 

k 

with the hamiltonian at sector k given as 

hk = -2tk3zk + 2Afcj2fc + 2efc (2) 

Here is the single electron energy, with = efc-f = e_fcj^, and is the pairing ("gap") 
energy . The dynamical symmetry , or spectrum generating algebra (SGA) , for each is 
su{2)k , so the SGA of the entire Hamiltonian is (g)^ su{2)k .When it causes no confusion we 
drop the index k. 

The su{2) pseudo-spin operators jpk obey 

[jpk^iqk]- = iepg^jrfc, where (p, g, r) = (1, 2, 3) (3) 

and each jpk is a bilinear in the fermion operators ( see Appendix 3). When all = 0, each 
hk is the hamiltonian for free electrons and the total wave function in the ground state is 
simply the product of individual creation operators acting on the vacuum state for that k. 
The filled Fermi Sea is the many-electron ground state , or the " disordered" state , denoted 
|0) = Da; o^tfcT^t-fcjO);.. For Afc 7^ the ground state is obtained at each k, by rotating hk 
about the jik axis by the angle 9ik = tan~'^{Ak/ek) , so that the transformed hamiltonian 
will be parallel to the j^k axis. 

The rotation operator is given by 

Uik = exp{i9ikjik) (4) 



so that 

hi = UikU{' = Ekjsk. (5) 
Now is directed along j^k and its eigenvalue is 



E, = ^Al + el (6) 

The key ingredient needed in order to discuss a phase transition is the order operator 
fjsc for superconductivity. A natural choice is the real part of the pair operator, i.e. fjsc = 
j2fc-Then, the order parameter in a state is the expectation value of that order operator in 
that state. We single out two states : the disordered state I'^dis) =|0), and the ordered state 
l^orcz)- Using the identifications given above for each of these states , depending on whether 
Afc = 0, or 7^ , we have rjsc = ,or rjsc 7^ respectively. We can go further . When 
Ak 7^ 0, the ordered state is l^iford) = Hkli^ord^k) = Ilfc ^ifc^'^lj"^) where \jm) is the ususal 
eigenstate of su{2) (for details see Appendix 3). Then, in the ordered state we have 

Vsc = (^fcb2fc|V'fc) (7) 

= {jm\Uij2kUi'\jm) (8) 

= {jml (j2fc cos 9i + jsk sin 9i)\jm) (9) 

= msin6'i (10) 

or 



r]sc = mA,/^Al + el (11) 

The lowest (ground state) energy occurs when m = —1/2, and the order parameter 
\Vsc\ 7^ in this state. 

The states \ipk) = U{^\jm), and in particular the ground state \ipkG) are su{2) coher- 
ent states in accord with the usual definitions [^], [0, ||3^. Globally the symmetry 



of the total hamiltonian H, which is ^ksu{2)k , gives a total wave function of the ground 
state which is the global coherent state 
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*BC5) = n IV'ikG) 
k 



(12) 



The steps used in this Section will be used again below. Namely, we: a) identify 
the dynamical symmetry of the hamiltonian chosen; b) diagonalize the hamiltonian ,thus 
obtaining the ground and excited state eigenfunctions (coherent states) and energies ; c) 
evaluate the expectation of the relevant order operator in the appropriate state to give the 
order parameter. 
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III. Algebra h(4) for a Displacive Ferroelectric Soft Mode. 

Ferroelectricity in perovskite-like systems is due to a "soft phonon" transverse optic 



(TO) lattice mode resulting from ionic displacements that break inversion symmetry pH 



||40|| , When these displacements are "frozen" by higher order anharmonic terms 
stabilizing a distorted structure, a macroscopic FE polarization P arises; the magnitude \P\ 
is proportional to the expectation value of the frozen soft mode amplitude (Qo)- Here Qo 
is the operator of the normal coordinate displacement for the soft mode In terms of 



the harmonic oscillator boson operators Bq, Bq of the soft mode, Qo ~ {Bq + Bq). We will 
assume that only a small region in the fc-space softens, and we can take that to be near 
k = . Thus we consider a homogeneous ferroelectric material , and we neglect FE density 
waves or stripes which may occur if A; 7^ 0. The soft mode will be represented by a harmonic 
oscillator. 

In the presence of spontaneous polarization P, a macroscopic self-electric field E 
arises , giving an energy proportional to E ■ P [^]. The field E is the internal transverse 
electric field generated by the transverse polarization . The origin of the latter is the frozen 
TO mode. Translating this self-term into the language of our problem we shall write this 
energy term as 7i£^(-Bq + Bq) , where 71 is a coupling constant we shall take as positive. 

We then take as our model for the ferroelectric sector of the hamiltonian the sum of 
the harmonic oscillator term for the soft mode , plus the energy due to the self-field coupled 
to the polarization ; higher order anharmonic terms are omitted. In second quantized form 
this Hamiltonian is : 

Hfe = iOToiBlBo + 1/2) + ^,£{Bl + Bq). (13) 

The frequency of the soft TO mode is uto, we take 71 as a positive constant , and £^ is the 
magnitude of E. We immediately recognize HpE as a" displaced oscillator" hamiltonian for 
the soft mode, including coupling to the self-macroscopic field. 

Before proceeding, note that the theory of ferroelectricity for soft mode perovskite-like 
systems has a long history. Microscopic models including anharmonic terms in the Hamil- 
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tonian followed the original work of Cochran (1959) [^] and Anderson who proposed 



the "soft mode" model . The models were carefully analysed by Cochran [Q, Cowley |[40|| , 
Cowley and Bruce [^], and others. Soon after experimental work on the doped tungsten 
bronzes was reported , several detailed theoretical papers appeared giving various micro- 
scopic models for the interplay between SC and FE , and also related to the Metal-Insulator 
transition in these materials |]l6l, [T^, |]l8l , and active study continues to the 



present 



Now we briefly comment about the soft mode. In the context of a traditional soft- 
mode displacive ferroelectric such as a perovskite like BaTiO^ or a tungsten bronze like 
A^WOs the soft mode can exhibit a typical temperature dependence such as: 

u;l, = n\iT-n)\ (14) 

with Q a constant , and Tc the FE transition temperature. On the other hand , in their 
very recent investigations on the mechanism of High Temperature Superconductivity in the 
cuprates, Weger and collaborators [|T|, 0, have examined quantitatively how the electron- 
phonon and electron-electron interactions are affected by the medium itself having a high 
lattice-induced dielectric coefficient. For example , values of e{uj) ~ 50 or more for u = 
lOmev , in La2-xSrxCu04 and YBCO were recently measured As noted by these 



authors, the relevant soft-mode transverse optic (TO) phonon is governed by the Lyddane- 



Sachs- Teller (LST) relation which for a single pair of LO and TO modes is: 



Here e(oo) is the high frequency dielectric coefficient, and ujlo,to are the characteristic 
phonon frequencies. This single mode expression works well for the c-axis component of the 
modes in both La2-xSrxCu04^ and YBCO . In our work we will not use the explicit expression 
given in equation (14) , but we shall make contact with the idea that the dielectric coefficient 
at low frequencies is very large , which is related to the "near ferroelectric instability". 

We return to the Hamiltonian for the ferroelectric sector, HpE It is well known |^ 
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|50| , that the "displaced oscillator hamiltonian" can be transformed by the unitary 



operator 

U2 = exp[UBl-Bo)]. (16) 

which displaces the oscillator Bose operators as 

U^Blu^' =Bl + io. (17) 

We take the simplest case with ,^0 real , and by choosing ^0 = (— 7i/t^ro) , "we obtain the 
transformed hamiltonian as 

H'pE = U2HfeU2^ (18) 

= [uto{bIb^ + 1/2) - (7i£)V^ro]. (19) 

This transformed oscillator hamiltonian Hp^ is shifted to a new minimum, but retains the 
same excitation frequency ujto as the original oscillator. 

We seek the eigenstates |$) of the original hamiltonian 

Hpe\^) = W\(^). (20) 

Transforming this equation by the operator U2 and continuing as in Section II , we find the 
eigenvalues and eigenf unctions as 

Wn = {n + l/2)-{^^S)yujTo. (21) 

and 

1$) = U-^\n) = exp - [UBl - Bo)]\n). (22) 

where |n) is a number eigenstate of the phonon number operator A''^ = BqBq. The state 
1$) is a Glauber coherent state for the FE oscillator [B? 



A natural choice of the order operator for the FE polarization is the coordinate 
operator Qq , or {Bq + Bq) . Thus fipE ~ <5o ~ {Bq + Bq). Clearly, in the state \n) i.e. , 
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I^dis) ) of the free phonon the expectation value of this order operator is zero: {n\Qo\n) — 
, meaning no spontaneous polarization . But in the Glauber coherent state |$) , which is 
|^ord)we have a non-zero value of the order parameter as the following argument shows: 

Vfe = i^m^) = {nm'iBi + Bo)U2\n) = 2^0 = -2{^iS)/u;to- (23) 

If, further, we identify the magnitude of the order parameter with the macroscopic 
polarization we have 

\Vfe\ = \P\ = {2-fi£)/ujTo. (24) 

Although the field E is a self-field due to the spontaneous polarization, we may treat this 
expression as defining a macroscopic dielectric susceptibility , given as x = {\P\/£) = 
2{ji/u!to). For a soft mode with ujto —>■ (for example as T ^ Tc) the susceptibility x will 
become very large . This interpretation agrees with what one would expect of a ferroelectric, 
or "near-ferroelectric" transition. Higher order terms will stabilize the system and prevent 
actual divergence. 

While the soft-mode displaced oscillator Hamiltonian is a very simphfied version of 
the true state of affairs, this hamiltonian HpE captures the physics of the FE sector for 
our purposes in this paper. Namely, this hamiltonian exhibits the FE displacement of the 
oscillator , which breaks a pre-existing inversion symmetry, and gives the enhanced dielectric 
susceptibility, in a simple algebraic h{4) setting. To include higher order anharmonic or 
coupled terms could make the model more realistic, but would depart from our algebraic 
framework. 

Just as we do not expect the su{2) SC model to be a microscopic model which 
can give all the features of superconductivity , so the /i(4) FE model does not claim to be a 
microscopic model incorporating interactions needed to explain all properties of ferroelectric 
media. 
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IV. Interaction Terms; Double Mean Field Approximation, 
a) Interaction Terms 

To precede we need the total Hamiltonian which we take as 

H = Hsc + HpE + Hjj^T- (25) 

We will take Hsc and Hfe as before, and now we turn to Hint ■ For the Hamiltonian 
to be translation invariant the interaction term must match the wave vectors of the 
soft-mode phonon and the Cooper pairs. Since the center of mass momentum of the pairs 
vanishes, we will directly couple the k = lattice soft mode and the electron pairs. 

In order to find the form of Hint in our algebraic framework, we need to use the basic 
operators in the SC and FE sectors, and combine them in an invariant fashion. It is natural 
to be guided by general prescriptions used in the Landau approach , and especially in 



the Ginzburg-Landau theory for competing order parameters. Here, for competing ferro- 
electricty and superconductivity, the order parameters are the spontaneous FE polarization 
\P\ and the superconducting gap A respectively. It is required , in the G-L theory ,that 
every term in the Free Energy shall be a scalar invariant under the relevant symmetry group, 
which in the present work is the direct product of the configuration space symmetry group 
of the crystalline medium, and the gauge group of the unbroken many-electron sector. We 
take the TO phonon in a prototype paraelectric crystal coupled to an "s-wave" SC complex 
gap parameter. For example take the prototype systems to be WO3 or perovskite-like , with 
cubic symmetry in the paraelectric and non-superconducting phase. The TO phonon will be 
split off from a r(15~) or other three-dimensional representation In the homogeneous 
approximation where the order parameters are uniform , the Free Energy density of the 
system will be of the form 

5F = a\P\^ + ib/A)\P\^ + a|A|2 + (3/A\A\^ + k/2\A\^\P\^ (26) 

In such a G-L theory , the lowest order, "generic" coupling between superconducting 
A and the ferroelectric P will be a biquadratic term of the form proportional to (|ApP^), 
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i.e. the last term in SF above. This is the lowest order term which satisfies gauge and parity 
symmetry requirements. A similar G-L theory was analysed by Liu and Fisher [Q, Imry 
||55| , Yurkevich, Rolov and Stanley |5^, and is discussed in texts such as that by Vonsovsky, 
Izyumov, and Kurmaev A rich phase diagram can result, depending on the relative 
magnitudes and signs of the coefficients of quadratic , quartic, and biquadratic terms .A 
case of particular interest to us, will occur for "antagonistic" order parameters for which k 
is positive. It is then possible to have four distinct stable regions in the phase diagram , 
including one for which there is coexistence of non-zero A and non zero P. Although our 
Hamiltonian will not contain terms of the fourth degree like \P\^ or |A|^, we will carry over 
the positive sign of the coupling constant denoted 72 below, corresponding to k. 

We need to translate this coupling term into the operators of our model. Thus we we 
must add to the free hamiltonian Hsc + , the interaction term , in which we use the 
correspondences |Ap ~ (j2)^ , and ~ Ql so that: 

Hint = 72(j1)(Q') = 72(j1)(^i + Bof. (27) 

We take 72, which corresponds to the K2, to be a positive coupling constant, as discussed 
above , in order to implement the " antagonism" of the two ordering fields. In principle both 
a microscopic theory of the basic interaction, and comparison with experiment ( see below 
in Section VI) will be needed in order to decide on the sign of 72 for a particular system. 

It should also be noted that if the initial paraelectric phase is of lower symmetry than 
cubic , a term of lower degree for example , |ApP could occur. This would preserve gauge 
invariance , and use the component of P transforming like the identity representation of 
the paraelectric group. Another manner in which a coupling linear in P might arise, uses 
the fact that the soft mode displacement, and the associated ferroelectric polarization P 
is transverse in the "long- wave" limit . So, in order to couple to P we should seek a 
transverse field operator associated to the superconductive sector. A natural candidate is 
the operator for the transverse superconducting current. Such a coupling term could have 
the form 
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^INT ~ ■ JsC (28) 

where the SC current operator would have the form 

J_5P = {-ie*h/2m)[^*V^ - ijVip*) - (e*Vmc) |VpA (29) 

where e* is the SC pair charge and A is the vector potentiaL In the London approximation 
Jsc ~ nsc-^ where nsc is the superconducting electron density , and V x A = h , with h the 
local magnetic field. Such a term which we can also write proportional to Pfe ■ A would be 
symmetry allowed in the correct geometry (lower than cubic , and with pre-existing broken 
i-symmetry). It would be linear in the transverse polarization. We could rewrite this term 
to bring it to the form of the operators of our model if we further take Jsc proportional to 
the number of pairs, as T,bl.bk and then express these in terms of the ji,2 of the SC sector, 
while P would be given in terms of the Bose creation and annihilation operators Bq, and 
Bq respectively. However , since these terms are linear in P and refer to pre-existent broken 
inversion symmetry , we do not include them here. 

The microscopic theory of the electron-soft ferroelectric mode (TO phonon) interac- 
tion was studied by Epifanov, Levanyuk and Levanyuk [^] . They point out that in the 



continuum approximation there is no transverse optic phonon to electron coupling because 
the TO mode does not induce a "macroscopic" electric field. By contrast, a longitudinal 
mode will give such a macroscopic field through the effective charge: V ■ P = —Ajipb where 



Pfe is an effective "bound" charge. But, as pointed out in |5^, using the correct lattice theory 



4^ , there will be coupling of the band electron to the local internal electric field due to the 



TO displacement. Taking the effective mass approximation , a meaningful continuum 
approximation can be , and was , defined for the electron-soft TO mode interaction and is 
discussed in Epifanov et.al. ||58|] . 

In very recent work Fay and Weger discuss the renormalization of the electron- 
phonon vertex in media with large dielectric constant near the ferroelectric phase transition. 

b) Double Mean Field Approximation 
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Assume now that the system is in a state such that the product of the squares of 
the fluctuations of the operators j2 , and Q evaluated in that state can be neglected. Thus 
if A and B are such operators and we write 

A=[A- {A)] + {A) = 5A + {A) (30) 

and similarly for then neglecting {{5AY) x {{SB)'^) , we will have: 

A^B^ ^ {2A{A) - {Af){2B{B) - {B^)) (31) 

Applying this to the biquadratic interaction term we obtain 

Hint - ]|(^o + ^o)' ~ "^APUbI + Bo) - 2AP% - 2PA\bI + Bo) + A^P^ (32) 

Our double mean field approximation (DMFA) yields a bilinear effective interaction term 
and it renormalizes the coefficients of the SC pairing term j2 and of the linear self-term 
in Q, and in addition there is an energy shift term proportional to |ApP^. In Appendix 1 
we examine the validity of the DMFA using the variationally determined wave functions. 
Recapitulating, we assume that initially we have the hamiltonian 

H^-2 ^{j^k + 2Akj2k) + u;o{Nb + 1/2) + ^iS{bI + Sq) + ^ l2k{j2k){Bl + Bo)\ (33) 

k k 

The bare coefficients A, and 71, refer to the prototype superconductor, and the prototype 
ferro-electric , and 72 is the initial pair-TO mode coupling coefficient. After making the 
double mean field approximation , and isolating a single mode k , we have the effective 
Hamiltonian at mode k in the DMFA : 

Hdmfa = -2ej3 + 2A'j2 + ojo{Nb + 1/2) + V^{bI + Bo) + V^UbI + Bo) + 72P'A2. 

(34) 

The renormalized coefficients are A', and ri 2. This will be our working Hamiltonian . 
It includes the SC and the FE prototype systems and their coupling via the soft-mode 
oscillator coupled to the pseudo-spin pairing Hamiltonian. Note that , technically , the 
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initial hamiltonian is in the "enveloping algebra" of su{2) ® h{4) because of the biquadratic 
terms , while Hdmfa is an element in the direct product algebra su(2) (8)/i(4). Extensions of 
our model to include various higher order terms coupling the phonon modes (anharmonicity) 
and pseudo-spin (pair-pair) terms will be treated elsewhere. 

The values of the renormalized coefficients in Hdmfa are , in terms of the original 
coefficients: 

A'^A(l-72n; 
Ti = 7i^-2(72PA2); 

T2 = 472. (35) 

When 72 ^ , we recover the sum of two separate sectors: for SC and FE. 

Our working model Hdmfa reduces to two other well-known models in certain limits: 
the " Jaynes-Cummings" model |61|, and the " spin-phonon" model |^2|, |^3[. The single 



mode "spin-boson" or "spin-phonon" model arises if A' = 0, and Fi = 0. Like our model, 
the spin-boson problem is not exactly soluble, but there has been much literature on it , 
including very recent work. Another related model is the single mode "Jaynes-Cummings" 
model for a photon (boson) coupled to a two- level "pseudo-spin" atom. We obtain the 
Jaynes -Cummings model by taking A' = , and Fi = , in Hdmfa , and if we write 

j2 =(l/2)(a+ + a_) = (l/2)(6t + 6) (36) 

where on the right hand side the and b refer to the electron pair operators, and , in the 
bilinear coupling term j2(-Bo +Bo) and if we retain only the "energy- conserving" interaction 
, we obtain: 

(a+Bo + cr^Bl). (37) 

This is the rotating-wave approximation (RWA) to the Jaynes-Cummings model. However, 
we cannot make this approximation for the SC-FE problem, that is for our HdmfAj as we 
will lose essential physics of our problem ; a similar point was made in the review of Leggett 
et.al. on the spin-phonon problem. 



V. The Variational "Coherent State" Eigenfunction. 

In order to obtain the ground state eigenfunction and eigenvalue of our model Hdmfa, 
we will use a variational procedure based on the su{2) (g) /i(4) symmetry. Recall that when 
72 = the uncoupled Hamiltonian at mode k is: 

H =hsc + HFE. (38) 

Recall that this Hamiltonian is " diagonalized" by the product of two unitary transformations 
previously denoted U — U1U2 where the parameters 6*1 = tan~^(A)/e , and ^0 = {—Ji/^^o) 
are the definite values obtained in Sections II and III . The resulting total state of the system 
without interaction |^) = U{^\jni)U2 ^\n) is a coherent state which is the product of the 
coherent state of the pseudo-spin su(2) algebra for the SC sector, times the Glauber coherent 
state for the h(4) algebra for the FE sector. 

For our coupled Hdmfa at mode k we introduce an analogous trial variational co- 
herent state (VCS) which is the product of two "coherent-like" states and is denoted |^^): 

1^,) = V\jm)\n) = V,-'\jm)V,-'\n) (39) 

where 

Vi = exp(i^ji) ;t>2= eMCiBl-Bo)) (40) 

but now, the parameters 9, and ^ are variational unknowns. The kets \ jm) and \n) are as in 
Sections II and III. We now define the energy in state as the diagonal value of Hdmfa 
in the variational coherent state |\E'„) 

(^.I^dmfaI*.) = E^,^{9,0 (41) 
and we determine 9 and ^ from 

dE^Jd9 = and dE^Jd^ = 0. (42) 

The evaluation of Emn is particularly simple since the only non-zero contributions to 
the diagonal value are from the matrix elements of the operators ja, Nb , and the constant. 
These enter as: 
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Csjs + ^qNb + Co 



(43) 



where 

C3= (2ecos^ + 2A'sine + 2r2Csin^). (44) 

and 

Co = uj^ie + 1/2) + 2V^i + TsC- (45) 

Hence, 

Emn{0,O^ mC3 + nojo + Co (46) 

and from this we have 

tane = -(A7e)-(er2)/e (47) 

and 

e = -Ti/coo - m(r2 sin 9)/u;o. (48) 

As a check of these results we verify that if 72 = then ^ — > ^1 , and ^ — > ^0 of Sections 2 
and 3 , with HipfT — 0- Corrections to the expressions for 9 and ^ are bihnear or quadratic 
in the couphng constants A', Fi , and T2, and are not considered here. 
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VI. Physical Consequences and Predictions, 
a) Energy and Order Parameters 

There are some immediate physical consequences of the previous results. First, con- 
sider the values of the SC and FE order parameters in the coherent state |^^) , where the 
variational parameters 6 and ^ take on the values just determined. 

For the SC order parameter we have 

= {3m\{n\V j2V~'^\3rn)\n) 

= (jm| (ja cos 9 - sin 9) \jm) 

= -msin6' (49) 

or 

r^lS' = ^(A'/e + er2/6)/[l + (A'/e + ^V,lefY'\ 
For the FE order parameter 

= {jm\{n\V{Bl + Bo)V-^\jm)\n) 
= (n|S^ + So + 2e|n) 
= 2^ 

or 

ril%^ = 2(-ri/cuo - r2(l - msin^))/cuo. (52) 

As a check we verified that in the absence of HiisfT,hoth order parameters rj^'^'^ revert to 
their values for the uncoupled systems with A' ^ A and ^ ^ ^o- 

Keeping terms up to second degree in A and P we find the order parameters in the 
coexisting phase . For the SC order parameter: 
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(50) 



(51) 



V^S' = Vscil - 12P') (53) 



A°(1-72P'). (54) 

7]FEil-l2^^)- (55) 
or 

In writing the above we identified the order parameters evaluated in the coexisting state 
with superscripts VCS, and we used superscipts for the gap, and polarization values with 
72 = 0, e.g. A° , and P° respectively, and we expressed the final results using the bare 
coupling parameter 72, instead of r2. 

At this point it is necessary to emphasize again that the sign of the parameter 72 
is not fixed by any symmetry argument, but must be determined from some microscopic 
considerations and comparison to experiments. As we pointed out in Section IV above, 
we have taken 72 positive in order to implement the competition between the two types 
of order. With this sign taken, we immediately conclude that: the presence of one non- 
zero order parameter will tend to supress the other one . Thus our model supports the 
"Matthias Conjecture" (1967) which we quote here: " Ferroelectricity seems to exclude 
Superconductivity more rigorously than Ferromagnetism seems to exclude Superconductiv- 
ity". We take the conjecture in a weaker sense for both coexisting superconductivity and 
ferromagnetism, and for superconductivity and ferroelectricity. Namely : in both cases there 
is an "adversarial tendency" , rather than some selection rule prohibiting coexistence. Since 
Matthias' statement, numerous examples of SC-FM and SC-AFM coexistence/competition 
have been found experimentally and studied theoretically. And, as pointed out in the Intro- 
duction, examples of the SC-FE coexistence are known in perovskite-type systems, and in 
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or 



^vcs 

For the FE order parameter: 



alkalie tungsten bronze doped systems. Older work on the P — W systems has considered 
superconductivity in "polar metals" ; some examples are in |12[, |T^, ||T^, |T^, p5| . 



Returning to our results , in the pseudo-spin model of pure superconductivity, the 
excitation energy at each k is given by E = [e| + A|]^/^. In our model , Hdmfa , the coupling 
will renormalize the gap parameter A to A'. So as a first approximation we can estimate 
the renormalized SC sector excitation energy and gap order parameter as 

A' = A(1-72P2). (57) 

or 

E' = E- (72 A2p2) ^ o(^2p4^ _ (5g^ 

The excitation energy E' is smaller in the presence of P 7^ 0, and the "SC gap", A' is also 
smaller than A. 

Similarly we note that the renormalized parameter Fi causes a shift of the mini- 
mum potential energy of the FE oscillator from —{'-fiS'^/uJo) to ~{'~fiS — 272PA^)^/ti;o, or 
~ — (7iP^/co'o + 47i72-EPA^/c(Jo. Thus there is a smaller downward shift of the potential 
energy of the FE oscillator, and a less stable minimum in the FE sector. 

Now consider the energy Emn{0,O which we have 

Emn = mCs + nuo + Cq. (59) 

We can identify the SC contribution as mCs, with m = (±1/2, 0) , which includes some FE 
admixture. Also (n + 1/2) — ^'^/ujq is the FE oscillator part, including some SC admixture. 
Such a separation is certainly not strict as the renormalized coupling constants A', Fi , and 
F2 involve all the interactions. 

b) Experimental Predictions 

Some experimental predictions follow from these results, 
i) Pressure Effect on T^. 
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From 

A^'='^ = A°(1-72P') (60) 

we can use the BCS result that 3.52kBT^ = 2A. We identify A° ~ T° as the "gap" 
or transition temperature when there is no FE. 

According to this resuh, the SC order parameter in the coexisting state, i.e. the 
"gap" A^'-^'^ , is decreased from A*^ by non- zero P . So, reducing P should increase the 
gap , and thus our prediction that: Tc should increase as P decreases. In order to test this 
prediction , we need a means of reducing P by some applied field. 

For most perovskite ferroelectrics, application of (positive) hydrostatic pressure vr will 
decrease P [41] . Thus dP/dn < 0,where vr is the applied hydrostatic pressure. If, further 
we assume that the effect of pressure on the bare gap (with no FE present) is small,i.e. 
dA^/diT ~ 0, then it is clear that for 72 > and dP/dn < we would have as one testable 
consequence of the above result 

dTjdn > (61) 

Experimentally, the pressure dependence of Tc in several Alkalie Tungsten Bronzes 
has been measured . Of particular interest for us are the sodium tungsten bronze family 
in which we distinguished the materials used by Matthias and others, with 0.1 < x < 1 , 
and the newer materials of Reich et.al. with x ~ 0.05. For the "Matthias type " sodium 
tungsten bronze : Nao^sWO^, it was found that 

dTc/dTT ~ +1.7 X W-^°K/bar (62) 

It is tempting to attribute this as confirmation of our prediction above. But no direct 
measurement of change of polarization with pressure vr was made. We encourage the mea- 
surement of the pressure- dependent ferroelectric polarization P(vr) in these tungsten bronze 
systems , which will enable a test of our prediction. 

We need also to recall other factors affecting the pressure dependence of the gap 



and Tc of superconductors. It is well known jS^] that a pure superconductor will exhibit 
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pressure dependence of Tc due to a number of factors : a) shift of the Fermi level under 
pressure, thus modifying the density of electron states at the Fermi level, b) change of phonon 
frequency under pressure ,and c) effects of pressure on defects , to mention some factors. 
(Additonal complexity is exhibited, for example, in V^Si, where hydrostatic pressure and 
uniaxial pressure in the [111] direction give positive coefficients, while uniaxial pressure in 
the [100] direction give a negative coefficient. [^ ). It is not simple to separate these effects, 
although for some specific cases theory was developed [Q, [|l^, [^, |6^. Also it is well 



known that even in the very well studied class of AB^ A-15, or (3 — W compounds , the sign 
of the pressure effect on Tc can be either positive or negative for different materials , with 
"no apparent universality" ||5^ . 



For the K and Rb tungsten bronzes the measured [Bloom et.al. ||6^ ] sign of the 
slope of pressure-dependence of is negative; however these materials have different crystal 
structures. As far as we can determine, there is no report of pressure effect on Tc for the 
materials used by Reich et. al., i.e. sodium tungsten bronzes with x ~ 0.05. 

To summarize : the available experiments on the older sodium tungsten bronze ma- 
terials of "Matthias type" [^ agree with our prediction . If the effect of pressure on the 



"pure" superconducting gap (no FE present) A° is small, then the experiments support the 
assumption that 72 > in our model hamiltonian. 

ii) Magnetic Field Enhancement of The Ferroelectric Polarization- A Non- 
Reciprocal Magneto-Electric Effect 

Turn now to the FE order parameter r]^^^ , which is the spontaneus FE Polarization. 

Using: 

P^^^ = P°(l-72|An, (63) 

where is the "bare" FE polarization in the absence of SC. For 72 > , P is decreased 
by the presence of |A| 7^ . Hence , the FE polarization should increase as A is decreased. 
Application of magnetic field will decrease A , ultimately to zero at He , the thermodynamic 
field. Hence we predict 
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dP/dH >0,H < He. (64) 

We are unaware of any experiments testing the predicted increase of the ferroelec- 
tric polarization P with applied magnetic field for the sodium tungsten bronze samples of 
Matthias type 0.1 < x < 1.0 used either by Matthias et.al. [1964, etc.] or by Bloom et.al. 
[1976] Nor have we found any reports of such measurements on other SC-FE materials, 
such as the doped SrTiOs, or the newer sodium tungsten bronzes. 

Our prediction of a change of the spontaneous polarization with applied magnetic field 
due to quenching of the superconductivity is a prediction of a new type of Magneto-Electric 
[ME] Effect, or more precisely, a Magneto- Polarization effect, since 

dP/dH ^ (65) 

signifies the Magneto-Polarization efffect in the coexistant SC-I-FE system. The usual 
Magneto-Electric Effect , predicted by Landau-Lifshitz [43], and first discussed for Cr203 
by Dzyaloshinsky, [^ is a property of materials which are usually ferromagnetic, or antifer- 
romagnetic, whose total symmetry group includes composite anti-unitary operations which 
will be broken such as rotation-reflection (or inversion) combined with time-reversal. Our 
coexistent SC-FE system breaks inversion plus gauge symmetry. Adding the applied mag- 
netic field to quench superconductivity breaks time reversal 6 , and places the SC-FE in 



presence of an external B field in a symmetry class for the ME effect [^], ||71|| . 

Thus a material in the coexistent superconducting-ferroelectric state ,in the presence 
of an applied external magnetic field will exhibit broken inversion plus time reversal and 
broken gauge symmetries. Hence our new Magneto-Polarization Effect is allowed. If we 
now use one of W. Pauli's famous aphorisms : " Anything not prohibited (by symmetry) 
is mandatory" we can then predict the existence on symmetry grounds of this new class of 
Magneto-Electric, or Magneto- Polarization Effect in a Superconducting Ferroelectric. There 
is, however one important remark to be made. Unlike the usual ME effect, in the SC-FE 
case, the effect is not "reciprocal" between magnetic field and applied electric field. 
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Previous work on the Magneto-Electric Effect has been reviewed authoritatively , 
||70| , [ fn]] . But, to our knowlege our prediction is the first of this new Magneto- Electric 
Effect in a Superconducting-Ferroelectric material. 
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VII . Discussion and Conclusions. 

The present work was motivated by our attempt to formulate the simplest algebraic 
model which would embody the relevant dynamical symmetries to describe the multicritical 
point for coupled superconductivity and ferroelectricity . In that framework it is natural to 
use the pseudo-spin su{2) algebra for superconductivity , since this algebra expresses the 
breaking of gauge symmetry embedded in the BCS theory [1957] [|^, and has proven useful 
in other contexts. We introduced the linearly displaced oscillator for the soft-mode phonon 
in describing a ferroelectric via the h{A) "Heisenberg Algebra" since this captures several 
key aspects of the displacive ferroelectric transition. Namely the shift to a new equilibrum, 
which breaks a pre-existing inversion symmetry and allows a spontaneous ferroelectric (and 
pyroelectric) moment, and the high dielectric coefficient. The simplest invariant coupling of 
these two order parameters which respects both the gauge and the inversion symmetries 
requires a biquadratic coupling which then turns to a bilinear coupling when reduced via 
our "Double Mean Field Approximation "[DMFA] . The Hamiltonian Hdmfa is an element 
in su{2) ® h{A). 

The physical idea motivating the present work is that the superconducting and fer- 
roelectric transitions are close to one another at a multicritical point. As we know from 
Ginzburg-Landau theory and renormalization group theories of such multicritical behavior, 
each transition will renormalize the other. A simple example of this is that our double mean 
field approximation (DMFA) renormalizes the "bare" coupling constants (see equation (35)), 
and leads to shifts in the two order parameters (see equations (60), (63)) from bare values. 
Hence the frustration, or reduction, of one order parameter will enhance the other. This is 
consistent with the weaker form of the "Matthias Conjecture" on the mutually antagonistic 
effects of superconductivity and ferroelectricity as shown earlier in this paper. Matthias was 
led to this statement in part by his work on the sodium tungsten bronzes. The new work 
by Reich et.al. on the sodium tungsten bronzes (in a different sodium composition range) 
reopens interest in the particular questions related to the superconducting-ferroelectric com- 
petition. We hope that our predictions in this paper and elsewhere on electrodynamics ||7^ 
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will stimulate further experimental work on these known coexistent systems including the 
titanate and bronze systems. 

Now we turn to possible relevance of this scenario to the high temperature cuprates. 
As was noted by Weger and collaborators [0]; conventional Eliashberg-McMillan theory 
of phonon-mediated interaction does not consider a multicritical point, nor in particular, 
the case where one phonon is a soft-mode ferroelectric phonon. There is no direct neutron 
scattering evidence for a soft ferroelectric phonon, however there is evidence that a close near- 
ferroelectric lattice instability exists in cuprates, and can play a role in the superconductivity 
mechanism. 

The measurements of ionic dielectric coefficients in several high cuprates shows 
exceptionally high values of the dielectric coefficients. Thus in YBCO and LaSrCuO val- 
ues of e ~ 40 — 50 were measured below dispersion frequencies of 19 mev and 27 mev, 
respectively. These values are reminiscent of the well-known near-ferroelectricity in the per- 
ovskites. [Recall that enhancement of SC by nearby FE was predicted for doped SrTiO^ 
in 1964 by M.L.Cohen |]73[, and was investigated thoroughly by A.Baratoff et.al.(1981) 



|7^. By itself doped SrTiO^ is not expected to be a superconductor, but rather a nearly 
ferroelectric semiconductor. However, the electron-phonon interaction is enhanced by the 
near-ferroelectricity, and possibly by the multivalley conduction band, and as a result doped 
SrTiOs is a superconductor with ~ 1.5°K.] 

In order to identify the phonons responsible for the high dielectric coefficient in the 
cuprates consider the dispersion of e{uj), which identifies the c-axis motion of the alkaline 
earth ion (Ba in YBCO, Sr in LaSrCuO) as responsible for the large values of e. These 
are not the phonons inducing the superconducting pairing , which are the planar oxygen 



displacements, at around 35 mev (transverse) according to |75[, and possibly near 70 mev 
(longitudinal) . 

We can then suggest the following picture which will relate to our model: the usual 
electron-acoustic phonon and possibly other electron-non-soft phonons are responsible for 



electron-pairing and their effect is present in the "pairing" coefficient A in Hsc 
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which carries over to the coupled Hamiltonian Hdmfa , which we use in our analysis. The 
separate electron-soft-mode or ferroelectric phonon represented by operators and Bq, 
with characteristic phonon frequency loto then couples to the formed electron pairs, giving 
the term Hint- 

In the language of Eliashberg-McMillan theory [^], |78| the electron-phonon coupling 
constant Aoo for coupling to phonon mode Qa is renormalized by proximity to a ferroelectric 
transition, so that the constant becomes frequency-dependent as X{uj) ~ e^(co')Aoo and 
also is large at low frequencies u -C ^Ia away from the physical sheet. This is consistent 



with the small phonon shifts in the superconducting state for , e.g. YBCO |^ where we 
can argue that X{uj) is very large. Note , however, large shifts were observed in BaCaHgCuO 
80| and in organic salts |81 . 



Making a short digression from the cuprates we can remark that the data on 
{BEDT — TT F)2Cu{C N 8)2 are illuminating in this context. The frequency of the phonon 
at 2 mev increases by almost 20 % below T^. The superconducting gap in this material 



is 2 A = lOmeV ||8^. According to the theory of Zeyher and Zwicknagl ||79|] , the phonon 
frequency should be pushed down below Tc, while experimentally it increases. We suggest 
that the mode at 2 mev is the soft or nearly "ferroelectric" phonon , which is not the one 
responsible for the pairing ; the phonons responsible for the pairing being around 6-8 mev. 
According to the results presented in this paper , superconductivity suppresses the tendency 
to ferroelectricity and thus , in effect , it "hardens" the relevant phonon. 

The data on HgBa2Ca^Cu40iQ is that the mode at 30 meV softens at by 
about 7 % , and the mode at 50 meV by about 10 %. Thus the direction of the effect is 



opposite to that in the organic superconductor. Hadjiev et.al. account for the effect as 
being due to the enormous anisotropy (which is the cause of the d-wave pairing). We do not 
consider this in the present work. Also , the frequency of the phonon is considerably higher 
than that of the "ferroelectric" mode , which is 19 mev in YBCO. 



The cited infra- Red measurements have shown 48 , [491, that the ionic dielectric 
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constant is dominated by the c-axis motion of the alkahne earth ion . Additional support 
for the existence of a large ionic dielectric coefficient associated with the motion of the Sr 
ion and its connection to the electronic properties of this material can be found from the 
EXAFS measurements of Pollinger et al ||8^ which show an anomaly in the distances of Sr 
to Oxygen in the LaSrCuO material. The large difference in distances between La-0 and 
Sr-0 supports a large polarization associated with the Sr. The interpretation is in terms of 
a Zhang-Rice singlet and an anti-Jahn Teller triplet which are nearly degenerate in energy. 
Some theoretical calculations of Anisimov -Andersen, and Kamimura , |^ support the 
near degeneracy ,which contributes to the high value of e. 

Another way to obtain the value of dielectric coefficient is via the small difference 
in energies between the Zhang- Rice singlet and the Jahn- Teller triplet. We can call this 
Ueff = Ubare/^ , which cau scrvc as a definition of e. Now, if we assume that this e is 
just that ionic dielectric coefficient actually measured in the IR experiments , then Ueff{uj) 
will be frequency dependent. Further, the small f/e// is associated with the lower edge of 
the "mid-IR" band , and it is described by a Jahn- Teller picture, as analysed by Moskvin 
et.al.(1998) |^. Added support for the identication of a small t/g// can come from the report 



that there is a softening of a longitudinal mode at a wave-vector near the zone boundary. 
Since this can indicate a tendency for charge segregation between the coppers it is consistent 
with a small U^ff- 

Returning to the sodium tungsten bronze systems there is a clear case for the ap- 
plicability of this model. Pure WO^ is ferroelectric , with large e, so the reported High 
gives a picture of a material which is a high Tc system in the presence of ferroelectricity 
but without copper which was previously anticipated to be ubiquitously associated with 
high temperature superconductivity. In the new bronzes now under study |p, the super- 
conductivity seems to be restricted to portions of the physical surface of the material ; this 
could be either an inherent effect, or possibly due to a preferential sodium concentration on 
the surface. The recent STM data of Levi, Millo et.al. indicates a superconducting gap 
which is sharp (as for s-wave pairing) but much smaller than in YBCO , suggesting weak or 
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intermediate coupling strength. Thus if the pairing is due to phonons , the frequency of the 
pairing phonon must be much higher-about 70 meV , as for example due to longitudinal 
oxygen displacements. Ferroelectricity here would involve displacements of the tungsten 
atoms, again consistent with two types of phonons playing different roles. 

In short, the cuprates may present a more complicated situation than the titanates 
or bronzes because the role of a latent or near-ferroelectric soft mode instability is partly 
masked - this is not surprising in view of the more complicated chemical composition and 
possible effects of magnetic excitations, as well as disorder in the cuprates. In the strontium 
titanate and sodium tungsten bronze cases application of the model seems straightforward, 
and it captures symmetry-related aspects of the coexistence-competition of the two collective 
effects. At a microscopic level our picture supports the view that there are two kinds of 
phonons playing a decisive role in these systems with SC-FE coexistence : first, the usual 
Frohlich coupled acoustic phonons giving the familiar BCS electron-electron pairing effects 
, and secondly, the , soft-mode ferroelectric phonons which yield the strongly enhanced 
static dielectric coefficient (e.g. via the Lyddane-Sachs- Teller related physics). As discussed 
elsewhere [Q, [Q, the latter phonons give an enhanced Thomas- Fermi screening length , 
which in turn changes the electron-phonon coupling and then the net effective electron- 
electron coupling, so that when inserted in the Eliashberg equations the new physics arises. 

The dynamical symmetry model proposed here captures certain essential features 
of the superconductor-ferroelectric competition/coexistence. Predictions of pressure and 
magnetic-field effects can be tested and will enable the sign of the coupling constant 72 
to be determined. For a given material the same coefficient will regulate the change in 
under pressure and the ferroelectric polarization under magnetic field as given in Section VI. 
Novel electrodynamic effects which are predicted for a nearly-ferroelectric superconductor 
are discussed elsewhere [0. 

VIII . Acknowledgements 

It is a pleasure to acknowledge helpful discussions and comments at various times 
from Professor Y. Avron , Professor C.W.Chu, Professor C.S.Ting, Dr. Mikhail Chernikov, 

32 



Professor M.E.Fisher, and Professor David Schmeltzer. Also help from Ms. Nguyen Que 
Huong is gratefully acknowledged (J.L.B.). J.L.B. acknowledges with thanks the hospital- 
ity of Profesor Ady Mann, and Professor Michael Revzen, and the Physics Department, 
Technion, Haifa, Israel, where part of this work was carried out. 



33 



Appendix 1: Calculation of Variances, and the DMFA. 

In order to examine the vaHdity of the DMFA we need to calculate variances of the 
operators j'l and Cf' = {Bq + Bq)^. First, we estimate the corrections to the parameters 
6 and ^ due to adding the interaction Hamiltonian Hint , which is characterized by the 
strength 72 , to the non-interacting hsc and Hfe- Going back to equation (54) and the 
renormalized coefficients given in equation (55) , after a little algebra , and retaining terms 
of lowest (linear ) order in 72 , we obtain: 

tan^ = (tan^i)(l - 72P^) (66) 

and , 

e = 6-4m72sin^i, (67) 

where sin 6*1 = {A/e)/^[l + (A/e)^]. Next, we note that j| = (l/4)(2fifcn_fe - 2j3) so that 
we will be able to easily find the expectation values directly. We then calculate 

VI ^ 2 cos^ OjI + (1/2) sin 29j2 + sin^ OjI (68) 

with Vi — exp{idji) as in equation (39). We then find 

{{5hr)^{{hY)-{{h)f^Scos'9 (69) 

and then as a figure of merit we can take: 

{iSj2r)/{{h)r = cot^^i(l + (l/2)72^P^sin^^0/(l - (70) 

Substituting the value of tan 9i , this normalized variance can then be expressed in terms of 
the parameters of the Hamiltonian as: 

Apart from the numerical factor , the scale of the variance is set by (A/e)^ which is clearly 
<^ 1. This will be true as well when the FE Polarization is non zero. Consequently we 
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verify that the variance of the operator representing SC ( the real part of the gap operator) 
is neghgable in the variational coherent state. 

Turning to the variance of Q , we need to calculate it using the operator {Bq + -Bq)^, 
but this is equal to {B^y + (Bq) + 2Nb + 1- Then the matrix elements are easily evaluated: 

{n\V2Q^Vf'\n) = 2n + 4(0' + 1. (72) 

and 

{n\V2QVf^\n) = 2^. (73) 

The normalized variance is then: 

HQ') - m?)im? = ((2n + i)/4(o'). (74) 

Since ^ ~ P/2 + 0(72) , and P will be large in the FE state , we verify that {{5Q'^)) ^ 1 in 
the variational coherent state |^^). 

We may then conclude that the DMFA , in the variational coherent state approxima- 
tion with wave function |^^) is a self-consistent approximation. 
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Appendix 2: Correlators 

Using our model we can examine the thermal average of two-time correlators of the 
basic operators of the model: the FE phonon displacement operator and the real part of 
the pair operator j2. The thermal average of the two time correlator of an operator Oa is: 

((C'a(i)C'a(O))) = {l/Z)Tr{exp - {(5H))exp{iHt)d^exp{-iHt)0^ (75) 

The partition function is Z = Tre~^^ with /? = {kBT)~^. Cross correlators can be similarly 
defined, by replacing one of the a subscipts by a different index referring to another of the 
operators. 

Calculation of these correlators is not possible if we use the full Hamiltonian, or even 
if we use Hdmfa- However if we denote the transformed Hdmfa by H' 

H' = VHdmfaV-^ (76) 

and evaluate H' at 6 — Oycs and ^ = ^vcs, we may then take as a good approximation 

i7' - Cs + cuoNb + Co (77) 

Making this Ansatz will permit us to easily evaluate the operators in Heisenberg 
picture,and then evaluate the traces. This is consistent with retaining the dominant leading 
order in the coupling parameters A, 71 , and 72. 

We first work out the displacement-displacement correlator. Using invariance of the 
trace we can then write 

((g(t)g(O))) = (l/Z)Tr[Ve-^^V-Ve*^*V-VgV-^ye-*^*V-VgV-^] (78) 

Using the properties of the operators Q and H transformed under V we obtain an 
intermediate result 

((g(t)g(O))) = (l/Z)rr[e*'^'^°^«(e*'"o*5(S + e-^^»*Eo + 2^) + 2^6"'^'^°^^ (g + 2^)]. (79) 

As only the lattice oscillator Bose operators have survived we now carry out the 
diagonal summation over the harmonic oscillator quantum number n. We then find 

36 



Z((3) = e^'""/^ X [2 smh(/3a;o)]~^ (80) 

and finally 

{{Q{t)Q{0))) = 4^' + e-^'^°* + 2Z(P) cosuot. (81) 

The same correlator can be evaluated for a Hamiltonian Hp^ — a;o(-^B+l/2) ,witliout 
the linear displaced term ; it gives the same result as above without the term 4^^. The value 
of C to 

be taken here is ^vcs from the solution of the variational problem. 
In exactly the same way , but now doing the trace over the 4 electronic states with 
m = 0, 0, ±1/2 we obtain: 

{{Ut)j2m) = (sin^ eves X cosh/3C3)[4(l + cosh/^Cs]"'- (82) 
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Appendix 3-Some su{2) results. 

In this Appendix we record some known results, which will be useful in various 
calculations throughout the paper. If (iki, are the electron annihilation/creation operators 

— * 

for wave vector k, spin (|), the relevent pair operators are defined by 

^ = 4t«-h; h = (biy-, fik = ait^feT- (83) 

— * 

The su{2) pseudospin algebra at each k is generated by 

jife = H/2)(6i-6,), (84) 

j,, = {l/2){bl + h), (85) 

j3ik = (-l/2)(nfc + n_fc-l) (86) 

so that 

dpkJqk]- = iep^rjVfe, where (p, r) = (1, 2, 3) (87) 

A basic set of states for each sector hk can be obtained by starting from su{2) eigen- 
states labelled by a pair of indices . Thus, since the eigenvalues of n±k are (0, 1) the states 
can be labelled by the eigenvalues (n^, n_jk). Or, we can use the (jm) labels of the ket |jm) 
referring to su{2) (suppressing k), which is an eigenstate of p and of j^. Here 

isljm) = m\jm); f\jm) = j(i + l)li"^)- (88) 

where we note that p — [nkfi^k + Jsfc + jffc)- Thus j'^ , as well as j'sfc , depends only on 
the number operators for ±/c and so their eigenvalues are easily computed from the allowed 
eigenvalues (nfe,n_fe). Also the same states are obtained by applying creation operators to 
the vacuum state like 0^10) etc. Enumerating these states at each k we have: 

(0,0)~|3/4,1/2)~|0) (89) 
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(0,l)~|0,0)~aL,jO) (90) 
(1,0) ~ ICQ) ~ 4^10) (91) 

(1, 1) ~ |3/4, -1/2) ~ al^a^jO) ~ 6l|0) (92) 

In each line, the states are presented in order, labelled by {nk,n_k) , by \j,m) , and 
finally by creation operators applied to the vacuum |0) . In constructing the eigenstates for 
hk these basic states of the free hamiltonian will be rotated to produce the su(2)-coherent 
states, which are eigenstates of hk- Recall that the states (0, 0) and (1, 1) are the basis for 



the irreducible representation D^^^'^\ and (0, 1) and (1,0) are bases for D^^^ , of s'u(2).(N.B. 
The latter two states are a time reverse pair for the corepresentation of su{2) (8> O). 

The eigenstates of j^k are the basic set \jm) given above. To obtain the eigenstates 
lipk) or I'^ord) , of the original problem 

kli^k) = Xkli^k) (93) 

transform hy Ui as follows 

UihkU^'Uil^k) = XkUim (94) 

and then we find 

Afc = mEk, and {ipk) = U^^\jm). (95) 
Now we can evaluate the order parameter in the state l^/'jt) = \^ord)'- 
rise = (V'fc|j2fc|V'fc) = {jm\Uij2kUi^\jm) = {jm\{j2k cos 9i + js^ sin6'i)|jm) = m sin 6^1 (96) 
or 



7^sc = mAk/jAl + el (97) 
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The lowest (ground state) energy Xq occurs when m — —1/2, and the corresponding order 
parameter \risc\ 7^ 0, and = {—l/2)Ek. The ground eigenstate \ipkG) is the BCS pair 
state which can be written in several alternate forms. These are 

l^kc) = exp{-ie,jik)\V^, -1/2) (98) 

^exp-[(e,/2)(bl-h)]bl\0) (99) 

= U-^\3/A,-l/2). (100) 

or 

\^kG) = -{^k + Vkblm (101) 

with Uk = —sin{0i/2) ; = cos{9i/2). 

In passing we note that we can recover the BCS gap equation from the above. The 
flsc given above is for one sector k. If we assume that = V/S. ,i.e. is independent of k 
and that the "global" order parameter is 

^ = Y.^scik) (102) 

k 

we may then obtain a self-consistent equation for the total gap function as: 



A = Y^^scik) = (1/2) Y.iV^,)/^Jel + (103) 

k k 

and substituting the above we have: 

2 = E(^)/V4 + A^ = V j p{e)de/V^^T^ (104) 

k 

Here,we changed variables to energy , and p(e) is the density of states (at the Fermi level) 
This reproduces the familiar result 

{2/V) = J p{e)de/VJT^. (105) 
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